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ABSTRACT 

We report the discovery of a brown dwarf that transits one member of the M+M binary system 
LHS6343AB every 12.71 days. The transits were discovered using photometric data from the Kelper 
public data release. The LHS 6343 stellar system was previously identified as a single high-proper- 
motion M dwarf. We use adaptive optics imaging to resolve the system into two low-mass stars 
with masses 0.370±0.009 M Q and 0.30±0.01 M Q , respectively, and a projected separation of 0'.'55. 
High-resolution spectroscopy shows that the more massive component undergoes Doppler variations 
consistent with Keplerian motion, with a period equal to the transit period and an amplitude consistent 
with a companion mass of Mq = 62.9 ± 2.3 Mj up . Based on our analysis of the transit light curve 
we estimate the radius of the companion to be Rq = 0.833 ± 0.021 Rj up , which is consistent with 
theoretical predictions of the radius of a > 1 Gyr brown dwarf. 
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1. INTRODUCTION 

Situated on the mass continuum between planets and 
hydrogen-burning stars are objects commonly known 
as brown dwarfs, with masses spanning approximately 
13Mj up up to 80 Mj U p (assuming Solar metallicity) . Since 
the first discoveries of these substellar objects fifteen 
years ago (Oppenheimer et al. 1995; Basri et al. 1996), 
various surveys have found additional examples in num- 
bers exceeding the population of known exoplanets 8 . 
However, despite the large sample of brown dwarfs, very 
little is known about the physical properties or forma- 
tion mechanism(s) of these substellar objects (e.g. Basri 
2006; Burgasser et al. 2007; Liu et al. 2008). 

Most known brown dwarfs have been discovered as soli- 
tary objects by wide-field, near-infrared (NIR) imaging 
surveys(e.g. Martin et al. 1997; Burgasser et al. 1999; 
Lawrence et al. 2007; Delorme et al. 2008). Their identi- 
fication is often based on spectral typing, with physical 
parameters derived from comparing photometric mea- 
surements to substellar evolutionary models. Knowledge 
beyond spectral typing is limited by the difficulty in mod- 
eling the complex molecular features that dominate the 
spectra of cool dwarfs (e.g. Allard et al. 2001; Cush- 
ing & Vacca 2006), a problem that persists above the 
hydrogen-burning mass limit for M-type dwarfs (Maness 
et al. 2007; Johnson & Apps 2009). 

The preferred method of measuring physical properties 
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of substellar objects such as masses, compositions and 
ages, is to study examples that are physically associated 
with brighter main-sequence stars. By assuming both 
the brown dwarf and its host star formed at the same 
time from the same molecular cloud, ages and chemical 
composition of the companion can be tied to the prop- 
erties of the brighter, more easily characterized compo- 
nent (Liu & Leggett 2005; Bowler et al. 2009). However, 
brown dwarfs in these favorable "benchmark" configu- 
rations are found in numbers far below the sample of 
exoplanets, despite their relative ease of detection com- 
pared to planet-mass companions. This observed feature 
of the substellar mass distribution of bound companions 
is known as the "brown dwarf desert," and the barren 
region extends over a wide swath around stars, extend- 
ing from w 0.05 AU out to hundreds of AU (Marcy & 
Butler 2000; McCarthy & Zuckerman 2004; Grether & 
Lineweaver 2006; Johnson 2009). 

The existence of a deep minimum in the mass con- 
tinuum between stars and planets suggests that distinct 
formation mechanisms operate at either mass extreme, 
one for stellar objects and one for planets. However, the 
scarcity of objects in the brown dwarf desert makes it 
difficult to determine where this line should be drawn. 
For example, it is unclear whether a 20 Mj up object in 
orbit around a main-sequence star formed like a mas- 
sive planet, or instead should be considered part of the 
extreme low-mass tail of the stellar initial mass func- 
tion (Kratter et al. 2010). Furthermore, an issue as 
fundamental as the mass-radius relationship below the 
hydrogen-burning limit is largely unconstrained by obser- 
vations. Understanding the nature and origins of brown 
dwarfs requires a much larger sample of detections. 

One promising avenue for increasing the sample of well- 
characterized substellar companions is through wide-field 
photometric transit surveys. Since the radii of objects 
are roughly constant from 1-100 Mj up (Baraffe et al. 
1998), transit surveys are uniformly sensitive to com- 
panions throughout the entire brown dwarf desert. The 
transit light curve, together with precise radial veloci- 
ties (RV), provide both the absolute mass (as opposed 
to minimum mass, M sinz) and radius of the companion, 
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thereby directly testing the predictions of interior struc- 
ture models (Burrows et al. 1997; Torres et al. 2008). 
Once transits have been discovered, the door is opened 
up a wealth of follow-up opportunities that can measure 
properties such as the brown dwarf's albedo, temper- 
ature distribution, emission spectrum and atmospheric 
composition (sec, e.g. the review by Charbonncau et al. 
2005). Further, studying the distribution of physical 
characteristics of companions, and their relationships to 
the characteristics of their host stars, can inform theo- 
ries of the origins of brown dwarfs in the same way that 
the statistical properties of exoplanets inform theories of 
planet formation (Fischer & Valenti 2005; Torres et al. 
2008; Johnson et al. 2010). 

In this contribution we present the discovery and char- 
acterization of a transiting brown dwarf orbiting a nearby 
low-mass star in the Kepler field. 

2. PHOTOMETRIC OBSERVATIONS 
2.1. Kepler Photometry 

The K epler space telescope is conducting a continuous 
photometric monitoring campaign of a target field near 
the constellations Cygnus and Lyra. A 0.95-m aperture 
Schmidt telescope feeds a mosaic CCD photometer with 
a 10° x 10° field of view (Koch et al. 2010; Borucki et al. 
2010). Data reduction and analysis is described in Jenk- 
ins et al. (2010b) and Jenkins et al. (2010c), and the 
photometric and astrometric data were made publicly 
available as part of the first-quarter (Q0-Q1) data re- 
lease. The Q0 data have a time baseline UT 2009 May 
5-11, and the Ql data span UT 2009 May 13 through 
2009 June 16. 

Among the 156,000 Long Cadence stellar targets in the 
Kepler field is the nearby, high-proper-motion M dwarf 
LHS6343 {a = 19' l 10 m 14.31 s , 8 = +46°57 m 25.0 s ; Reid 
et al. 2004). The photometric properties as listed in 
the Kepler Input Catalog (KIC; Batalha et al. 2010) are 
given in Table 1. The Q0-Q1 photometric time series of 
LHS6343 contains a total of 2115 brightness measure- 
ments with a 29.4 minute cadence and a median internal 
measurement precision of 7 x 10~ 5 (Jenkins et al. 2010c). 

As part of a study of the photometric variability of the 
closest stars in the Kepler field, one of us (K.A.) noted 
that the light curve of LHS 6343 exhibits four deep, pe- 
riodic dimming events spaced by 12.71 days. The light 
curve depths are constant to 0.3 mmag and exhibit no ob- 
vious additional dimming at intermediate periods, con- 
sistent with the signal of a transiting planet-sized ob- 
ject. The astrometry shows no shift in the center of light 
greater than 1 millipixcl (4 mas, Jenkins et al. (2010a)), 
and there are no secondary eclipses evident at intermedi- 
ate phases. It is therefore unlikely that the source of the 
dimming events is a background eclipsing binary (EB). 

Examination of the Palomar Observatory Sky Survey 
(POSS) images shows no other stars at the current posi- 
tion of LHS 6343, further ruling false-positives involving 
an EB. The closest star in the archival images is 7'.'0 to 
the West of LHS 6343, near the edge of the Kepler photo- 
metric aperture. However, the star is only 4% as bright 
as LHS 6343, meaning that if it is an eclipsing binary it 
would have to nearly disappear to replicate the observed 
transit signal. This situation is ruled out by the lack of a 
large photocenter shift seen in the astrometric measure- 



TABLE 1 

Observed Properties of LHS 6343 



Parameter 


Value 


Source 


(X 


19 10 14.33 


KIC 
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+46 57 25.5 


KIC 


H-a. yi j 


-145 


KIC 


fig (mas yr ) 




1 \ IV. 


g' 


14.03 ± 0.02 


KIC 


r' 


13.06 ± 0.02 


KIC 


■t 
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i o Q7 _l n no 


IV 1L 


Btot 


15.009 ± 0.025 




Hot 


13.435 ± 0.018 




-Kp.tot 


1 Q 1 C\A _l_ f\ C\A 

13.104 ± 0.04 




T 

Jtot 


9.5 IV) ± 0.021 


oi\ IAOC 




8.972 ± 0.027 


2MASS 


Ks,tot 


8.695 ± 0.011 


2MASS 


AJ 


0.49 ± 0.05 


PHARO 


AH 


0.48 ± 0.05 


PHARO 


AK S 


0.45 ±0.06 


PHARO 


Ja 


10.10 ±0.04 




Jb 


10.59 ±0.06 




Ha 


9.51 ± 0.04 




H b 


9.99 ± 0.07 




K S ,A 


9.25 ± 0.05 




K S ,B 


9.70 ± 0.08 





ments. The Kepler photometric measurements phased at 
the 12.71 day period are shown in Fig. 1. 

2.2. Nickel Z-band Photometry 

We observed the transit event predicted to occur on 
UT 2010 June 29 using the 1-m Nickel telescope at Lick 
Observatory on Mt. Hamilton, California. We used 
the Nickel Direct Imaging Camera, which uses a thinned 
Loral 2048 2 -pixel CCD with a 6.3' square field of view 
(Johnson et al. 2008). We observed through a Gunn 
Z filter, used 2x2 binning for an effective pixel scale 
of C/37 pixel -1 , and a constant exposure time of 75 sec- 
onds. We used the slow readout mode, with 34 s between 
exposures to read the full frame and reset the detector. 
The conditions were clear with <~ l'.'O seeing. We began 
observing as soon as possible after sunset at an airmass 
of 1.4 and observed continuously for 5.4 hours bracketing 
the predicted transit midpoint, ending at an airmass of 
1.05. 

We measured the instrumental magnitude of LHS 6343 
with respect to the four brightest comparison stars in 
the field using an aperture width of 23 pixels and a sky 
annulus with an inner and outer radius of 28 and 33 pix- 
els, respectively. We converted the Nickel timestamps to 
BJDjjtc using the techniques of Eastman et al. (2010) 
to be consistent with the Kepler data. The Nickel pho- 
tometric measurements phased at the 12.71 day period 
are shown in Fig. 1. 

3. RADIAL VELOCITIES AND ORBIT SOLUTION 

We obtained spectroscopic observations of LHS 6343 at 
Keck Observatory using the HIRES spectrometer with a 
resolution of R ~ 55, 000, with the standard iodine-cell 
setup used by the California Planet Survey (Howard et al. 
2010). The transit depth, together with a rough stel- 
lar radius estimate of 0.4 R Q appeared consistent with 
a planet with a radius of ~ 0.7 R,j up . In anticipation 
of a low-amplitude Doppler signal we initially used 45 
minute exposures through the iodine cell and C2 decker 
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Fig. 1. — The Kepler (upper) and Nickel (blue) light curves, 
phased at the photometric period. The best-fitting light curve 
models are shown for each data set (see § 5.3), and the residuals 
are shown beneath each light curve. The Kepler residuals have 
been multiplied x5 for clarity. 



for sky-subtraction. The resulting signal-to-noise ratio 
(S/N) was s» 90 at 5500 A, near the center of the iodine 
absorption region. 

A cross-correlation analysis of the first two observa- 
tions revealed two peaks separated by ~ 10 km s" 1 . 
The lack of a coincident background star in the POSS 
images suggests that the second set of lines must be 
from a physically-associated binary companion and that 
LHS 6343 is a double-lined spectroscopic binary. Adap- 
tive optics observations described in § 4 confirmed the 
existence of a wide binary companion at a projected sep- 
aration of 0'.'55. Hereafter, we refer to the more massive 
component as Star A, and the less massive component 
as Star B. 

In order to discern which component of the binary sys- 
tem is transited by a companion we made subsequent 
HIRES observations with a position angle oriented along 
the binary axis to ensure the light from both stars fell 
within the slit. The cross-correlation analysis of the third 
observed spectrum revealed that the deeper set of ab- 
sorption lines shifted by ~ 10 km s _1 with respect to the 
first observation, indicating that Star A is orbited by a 
massive companion. 

Our remaining HIRES spectra were obtained without 
the iodine cell and with 3-minute exposure times. We 
measured the radial velocity of Star A using the cross- 
correlation analysis described by Johnson et al. (2004). 
However, we modified the procedure by constructing a 
double-lined cross-correlation template. We began with 
an iodine- free spectrum of the HD 265866 (M3V) for star 
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Fig. 2. — Keck/HIRES radial velocity measurements of 
LHS 6343 A. The best-fitting Keplerian orbit solution is shown as 
a dashed line. The systemic velocity, 7 = —46.0 ± 0.2 km s _1 , has 
been subtracted for clarity. The lower panel shows the RV residuals 
about the best-fitting orbit. 



A, and added to this template a scaled, shifted version 
of itself to represent the spectrum of star B. For each 
observation we adjusted the scaling and Doppler-shift of 
spectrum B with respect to spectrum A until the cross- 
correlation peak was maximized. We then measured the 
centroid of the optimized cross-correlation function by 
fitting a parabola to the region near the resulting single 
peak. 

We corrected for shifts in the HIRES detector by using 
the telluric lines in the 630 nm a-band as a wavelength 
reference. We measured the absolute radial velocities of 
LHS 6343 A with respect to the absolute radial velocity of 
HD 265866 (M3V; V r = +22.97 ±0.50 km s" 1 ; Chuback 
et al., in prep), and the resulting RV is corrected to the 
Solar System barycenter using the velocity corrections 
computed by the CPS data reduction pipeline. The full 
time series of radial velocity measurements is displayed 
in Fig. 2. The radial velocities are also listed in Table 2, 
along with the heliocentric Julian Dates of observation 
and internal measurement errors. 

We searched for the best-fitting Keplerian orbit solu- 
tion using the partially linearized, least-squares fitting 
procedure described in Wright & Howard (2009) and im- 
plemented in the IDL software package RVLIN 9 . We fixed 
the period and mid-transit time based on the light-curve 
analysis in § 5.3, which leaves only four free parameters: 
the velocity semiamplitude (Ka), argument of perias- 
tron (w), systemic velocity (7) and eccentricity. We find 
that the RVs are described well by a nearly circular or- 
bit (e = 0.056 ± 0.032) with a velocity semiamplitude 
K = 9.6 ± 0.3 km s _1 . The full spectroscopic orbit is 
given in Table 5 and shown in Fig. 2. 

The parameter uncertainties were estimated using a 
bootstrap Monte Carlo algorithm. For each of 5000 real- 
izations of the data, the measured RVs are perturbed by 
adding residuals randomly drawn from about the best- 
fitting orbital solution, with replacement. We chose this 
technique over an MCMC analysis out of concern that 
our RV measurement uncertainty is dominated by sys- 

9 http://exoplanets.org/codc/ 
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Fig. 3.— PHARO if s -band adaptive optics image of LHS6343, 
showing the two M-type components of the system. Star A is to 
the upper right of the image and Star B is to the lower left. 

tematic errors related to imperfect treatment of the sec- 
ond set of absorption lines, rather than photon noise. 
Thus, instead of assuming the RVs are normally dis- 
tributed about the model, we use the residuals them- 
selves as an estimate of the noise model. 

4. PALOMAR ADAPTIVE OPTICS IMAGING 

We acquired near-infrared images of LHS 6343 on UT 
2010 June 29 using the Palomar 200-inch telescope adap- 
tive optics system and PHARO camera (Hayward et al. 
2001; Troy et al. 2000). These diffraction- limited ob- 
servations spatially resolve the target into a binary, as 
shown in Fig. 3. We used these observations to calculate 
the relative brightness of each component of the visual bi- 
nary and to subsequently constrain the range of possible 
masses (see § 5). Our results for differential magnitudes 
in each of the J,H,Ks filters are listed in Table 1. 

We used aperture photometry To measure the differen- 
tial magnitudes of the two stars. However, given their rel- 
atively small angular separation, special care was taken 
to account for cross-contamination between the compo- 
nents. To remove the majority of flux contributed by the 
neighboring star, we used the spatial symmetry in the 
images. The amount of contaminating starlight was esti- 
mated by summing the counts over a region the same size 
as the photometric aperture located on the side opposite 
the star of interest. 

Once the contaminating light from the neighbor star is 
removed, we find that the photometric precision is lim- 
ited to several percent by uncertainties resulting from the 
subtraction residuals, as well as CCD non-linearity and 
PSF centroiding. These errors sources contribute simi- 
larly to the overall uncertainty and were added in quadra- 
ture, neglecting any correlations between PSF centroid- 
ing errors and contamination removal, which we found 
to be comparatively small. To this, we also added in 
quadrature the standard deviation in the mean flux ra- 
tio of the companions over the 20 images acquired in each 



TABLE 2 
Radial Velocities for LHS 6343 



JD 


RV 


Uncertainty 


-2440000 


(km s -1 ) 


(km s- 1 ) 


loo/o.Uyo 


-38.40 


0.21 


15373.998 


-37.51 


0.17 


153/ I.VIo 


-49.42 


0.50 


15377.099 


-49.56 


0.45 


15378.030 


-52.56 


0.41 


15379.052 


-56.39 


0.44 


15380.827 


-55.83 


0.48 


15380.831 


-54.84 


0.47 


15395.983 


-47.92 


0.57 


15396.970 


-42.63 


0.51 


15404.974 


-55.90 


0.59 


15405.821 


-56.32 


0.46 


15406.865 


-53.69 


0.51 


15407.854 


-51.31 


0.57 



bandpass. Observations in the Ks band have a slightly 
larger uncertainty than in J and H, since the binary sep- 
aration subtends a smaller angle on the sky in units of 
resolution elements. 

5. STELLAR PROPERTIES 

The physical properties of the two stellar components 
of the wide binary, hereafter Star A and Star B, are of 
central importance to measuring the properties of the 
substellar companion, LHS 6343 C. The mass of the com- 
panion is related to the mass of Star A (Ma) and the 
companion's radius depends on the stellar radius, Ra- 
The luminosity of Star B is also important for the mea- 
surement of the companion's radius, as its contribution 
to the total flux of the system dilutes the transit depth. 
In fact, for the specific case of LHS 6343 the precision 
with which we can measure the companion radius will 
depend critically on our estimate of the "third light" con- 
tribution of Star B, rather than the photometric precision 
of the transit light curve (Irwin et al. 2010). 

Because low-mass stars spend nearly their entire lives 
close to the zero-age main sequence, their observed prop- 
erties are a function of two physical characteristics: mass 
and, to a lesser extent, metallicity. There are two widely- 
used methods estimating the masses of M dwarfs. The 
first involves a comparison of observed properties such 
as absolute magnitude and color index, or luminosity 
and effective temperature, to tabulated stellar evolu- 
tion models, such as those calculated by Baraffe et al. 
(1998). However, studies of low-mass eclipsing bina- 
ries (EB) have demonstrated that these models system- 
atically under-predict stellar radii (Ribas 2006; Lopez- 
Morales 2007), even in cases when stellar activity should 
play a minimal role in shaping stellar structure (Torres 
2007). 

The other method makes use of empirical relationship 
between the near-infrared luminosity of a star and its 
mass, as parametrized by Delfosse et al. (2000). We use 
the empirical relationships almost exclusively in order to 
avoid any systematic errors in the stellar evolution mod- 
els. However, since the mass-luminosity relationships re- 
quire absolute magnitudes, the distance to the star must 
be known to accurately estimate the stellar mass. Unfor- 
tunately, there is no published trigonometric parallax for 
LHS 6343, and the spectroscopic parallax of Rcid et al. 



5 



TABLE 3 

Polynomial Coefficients for Magnitude-Mass 
and Mass- Radius Relationships 



3 


bj 


b H 


b K s 


bv 


b R 





14.888 


13.211 


13.454 


-9.229 


0.000 


1 


-74.375 


-57.464 


-67.439 


5.017 


1.268 


2 


376.72 


271.11 


338.43 


-0.6609 


-1.013 


3 


-1089.6 


-762.12 


-976.27 


0.03314 


0.9391 


4 


1601.6 


1103.9 


1433.9 






5 


-935.37 


-640.80 


-838.03 







(2004) is unreliable because it is based on the total mag- 
nitude and colors of the binary system, rather than the 
individual stars. 

While the binarity of LHS 6343 in some regards poses 
a challenge, having two stars with the same age and 
chemical composition, together with the available pho- 
tometric measurements, provide a unique opportunity to 
determine the physical characteristics of the two compo- 
nents. As we will demonstrate, the luminosity difference 
between the two stars constrains the mass ratio, while 
the total luminosity constrains the total mass, distance 
and mctallicity of the system. 

An additional constraint on the mass is provided by 
the shape of the transit light curve. The slope of the 
ingress/egress yields the scaled semimajor axis cl/Ra = 
a_R, which is related to the density of Star A through 
Kepler's third law (e.g. Seager & Mallen-Ornelas 2003; 
Sozzetti et al. 2007; Winn 2008). However, the true 
value of an depends on the true transit depth (Rc/Ra) 2 , 
which in turn is related to the amount of dilution in the 
light curve due to Star B. 

In what follows we first relate the masses of Star A and 
Star B to the available observables. We then estimate 
the flux contribution of Star B in the Kepler bandpass, 
which provides a refined estimate of the stellar masses. 
The iteration of this procedure yields the stellar param- 
eters of Star A, which allows us to estimate the physical 
characteristics of the transiting object LHS 6343 C. 

5.1. Stellar Masses 

The most useful data available to us comprises seven 
photometric measurements: the total near-infrared mag- 
nitudes in the 2MASS catalog denoted by % where i = 
J, H, Ks corresponds to the three bands, respectively; 
the magnitude differences A J, AH, AK$ from our AO 
imaging, denoted by Aj. We also have the total John- 
son V- and £>-band magnitudes V to t — 13.435 ± 0.018 
and B tot = 15.009 ± 0.025. Our V tot agrees well with 
the value listed in the TASS catalog, V tot = 13.38 ± 0.24 
(Droege et al. 2006). These magnitudes yield the system 
color (B - V)tot = 1-574 ± 0.031. 

The individual apparent magnitudes in the ith NIR 
band (JHK) of Star A (to^a) and Star B (m^s) are 
related to the total magnitudes % and magnitude differ- 
ences Aj through 



™ M = 2.51og 10 (l + 10°- 4A ')+7; 
TOjR = A, + mi, A 



(1) 



Eqns 1 give three NIR apparent magnitudes for each star. 
These six NIR photometric measurements can be related 



to the stellar masses (Ma, Mb) through the following 
equations 



mj(M, d)=Mi(M) + 5(log 10 d - 1) 



(2) 



The functions Mi(M) give the absolute magnitudes in 
the NIR bands as a function of stellar mass, M, as de- 
termined by inversion of the empirical relationships of 
Delfosse et al. (2000), which we approximate with the 
polynomial 



3=0 



(3) 



The coefficients {bi} are listed in Table 3 for the i = 
{J, H, K s } bands. 

Under the assumption that the binary components 
share the same chemical composition, they should reside 
at the same distance from the average main sequence 
in the {V - K Sl M Ks } plane (Johnson & Apps 2009, 
hereafter JA09). JA09 provide a relationship between a 
star's metallicity, [Fe/H] = F, and its "height" above the 
Solar-neighborhood mean main sequence, AMk- Since 
the stars share the same composition the must lie on the 
same isometallicity contour: their V — Ks colors must 
be consistent with the same value of AMk, while the 
individual U-band luminosities must be consistent with 
the measured Vtot- This constraint can be expressed as 

V tot (M A ,M B ,d,F)=M Y (M A , F) 

-2.51og 10 (l + iqO-Hm v (m a ,f)^m v (m b ,f)]^ 
+ 51og 10 d-5, (4) 

where the absolute V-band magnitude, Aiy(M, F), is 
related to stellar mass M and metallicity F by inverting 
the photometric metallicity calibration of JA09 10 . We 
approximate this inversion using the polynomial 



M V (M, F) 



3=0 



Mks(M) 



F-0.05 
0.55 



1 1 



(5) 



The coefficients {by} arc listed in Table 3. 

In addition to the apparent magnitudes, the tran- 
sit light curve provides an additional constraint on Ma 
through the scaled semimajor axis cl/Ra = ifi- This 
quantity is related to the mass and radius of Star A, and 
the period and mass of LHS 6343 C, through 



a R (M A ,M c ,P). 



G 
4^2 



1/3 



M 



1/3 



Ra(Ma)' 



.p2/3 1 



M c 
M A 



1/3 



(6) 



10 We use the JA09 relationship rather than the Schlaufman & 
Laughlin (2010, SL10) calibration because the former provides a 
better match to the mean metallicity of the Solar neighborhood. 
In principle, the SL10 relationship would serve our purposes just as 
well since both V-band metallicity relationships provide a means 
of relating Vtot to the stellar masses under the constraint that both 
stars lie on the same isometallicity contour in the {V — Kg, Mjf s } 
plane. The only difference is the exact value of [Fe/H], which will 
be ~ 0.1 dex lower using the SL10 calibration. 
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where R A (M) is a function relating the stellar radius and 
mass. We use a polynomial fit to the masses and radii of 
well-characterized {{(Jm,&r\ < 3%) low-mass eclipsing 
binaries tabulated by Ribas (2006) 



Ra(M) =J2bR,jMi 



(7) 



1=0 



The coefficients {bp} are listed in Table 3. 

The optimal stellar parameters can be obtained by 
minimizing the fitting statistic 



3 2 

2 _ ( ™>i,A - rrii(M A ,d) 

Xtot — 



1=1 

3 

+E 

»=i 



rrn.B - mi(M B ,d) 



+ 



+ 



V tot -V tot (M A ,M B ,d,F) 



°~V tc 



a R -a R (M A ,M c ,P) 



(8) 



However, in order to obtain the best-fitting parameters 
we must determine the corrected values of the transit 
depth and ap by accounting for the flux contribution of 
Star B, as described in the following section. 

5.2. Flux Contribution of Star B 

Since both stars fall within the Kepler and Nickel pho- 
tometric aperture, we must estimate the flux contribu- 
tion of Star B in both the Kp and Z bandpasses. The 
transformation between the Johnson B and V magni- 
tudes and the Kepler magnitude Kp is given in the Ke- 
pler Guest Observer web page 11 . However, while we were 
able to measure individual V magnitudes for both stars 
based on their total V magnitude and derived physical 
properties using Eqn. 5. However, there is no suitable 
relationship between the available observables and the 
individual B magnitudes. We are therefore forced rely 
on stellar model grids to estimate the B — V colors of 
stars A and B. 

For this task, we selected the "Basic Set" of stellar 
model grids from the Padova group for Solar composition 
and an age of 5 Gyr (Girardi et al. 2002), which give 
predictions for the Johnson B and V magnitudes as a 
function of stellar mass. We selected the Padova models 
because the Baraffc et al. (f 998) model grids do not give 
fluxes in the B or Z bandpasses. We tested the reliability 
of the Padova model grids using a sample of 6 metal-rich 
stars from Johnson & Apps (2009). We found that the 
models under-predict the B — V colors of our calibration 
stars by 0.20 ± 0.07, independent of stellar mass over a 
range of approximately O.f M Q to 0.6 M©. We applied 
this correction to the Padova models to obtain a refined 
relationship between stellar mass and B — V color, and 
the uncertainty of this offset was folded into the final 
uncertainties in the colors. 

11 http://keplergo. arc .nasa.gov/CalibrationZeropoint . shtml 



Using the estimates of the B ~ V colors of Star A and 
Star B, together with the individual U-band magnitudes 
from Eqn. 5 we can estimate the relative magnitudes of 
the stars in the Kepler bandpass, AKp, as a function 
of the stellar masses. We then use this value of AKp to 
correct the Kepler light curve parameters. For the Nickel 
light curve, we use the SDSS z-band magnitude from the 
Padova model as a proxy for the Gunn Z filter. 

Our reliance on low-mass stellar model grids for the 
magnitude differences is less than ideal, particularly in 
optical bandpasses, because of the difficulty in properly 
treating molecular opacities at optical wavelengths (e.g. 
Baraffc et al. 1998). In the analysis that follows, we at- 
tempt to ameliorate this imperfect knowledge by using 
large uncertainties for AKp and AZ, which are propa- 
gated throughout our analysis and reflected in the con- 
fidence intervals for our derived system properties. We 
discuss the impact of our model-based magnitude differ- 
ences in § 6.2. 

5.3. Light Curve Analysis 

We fitted the Kepler and Nickel light curves using the 
analytic eclipse model of Mandel & Agol (2002) to com- 
pute the integrated flux from the uneclipsed stellar sur- 
face as a function of the relative positions of the star 
and planet. The parameters were the period P, inclina- 
tion i, ratio of the companion and stellar radii Rc/R A , 
the scaled semimajor axis ap, time of mid-transit Tmid) 
and the parameters describing the limb-darkening of the 
star. For the limb-darkening we used the quadratic ap- 
proximations of Sing (2010) and Claret (2004) for the 
Kepler and Z bands, respectively. In our model we fixed 
the quadratic terms u 2 at the tabulated values for each 
filter, and allowed the linear term u\ to vary under a 
penalty of the form exp[— (x — /x^)/<r^], which is added 
to the fitting statistic. This treatment of the limb dark- 
ening was made based on the additional structure seen 
in the residuals of the Kepler light curve fit. Allowing 
the limb-darkening parameters to float freely results in 
unphysical values. 

To properly model the light curves we made two mod- 
ifications to the typical light curve analysis. First, we 
corrected for the third-light component by adjusting the 
normalized flux level from the Mandel & Agol light curve 
model, fmod, such that 



fcA 



mod 



1 + IQ-OAAKp 



(9) 



In our fitting procedure, we treat AKp as a free pa- 
rameter under the normally-distributed, prior constraint 
JV(AKp,aK P ) based on the estimate of the magnitude 
differences described in § 5.2. We use a similar prescrip- 
tion for the flux contribution of Star B to the Nickel Z 
bandpass. Our other modification involved rebinning the 
analytic light curve to match the 29.4-minute cadence 
of the data using a method similar to that of Kipping 
(2010). 

In addition to modeling the companion transit, we also 
fitted a slowly- varying function to the out-of-transit por- 
tion of the Nickel light curve to account for differential 
extinction from the Earth's atmosphere. To account for 
time-correlated noise in the data we used a Daubechies 
fourth-order wavelet decomposition likelihood function 
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TABLE 4 
LHS 6343 Ephemeris 



T mid (BJD - 2450000.0) 


T m id - Ephemeris 


Telescope 


4957.216636 ± 0.000073 


0.0000020 ± 0.00012 


K 


4969.930410 ± 0.000099 


-0.000045 ± 0.00014 


K 


4982.64437 ± 0.00013 


0.000093 ± 0.00016 


K 


4995.35807 ± 0.00015 


-0.000028 ± 0.00018 


K 


5376.77274 ± 0.00025 


-0.0000050 ± 0.00027 


N 



Note. — K - Kepler. N - Nickel Z-band 

following the technique described by Carter & Winn 
(2009). Wavelet decomposition provides increased confi- 
dence in the derived parameter uncertainties over the tra- 
ditional x 2 likelihood by allowing parameters that mea- 
sure photometric scatter (uncorrelated white noise a w , 
and 1// red noise ay) to evolve as free parameters. The 
method recovers the traditional \ 2 fitting statistic in the 
case of ay =0, and when a w is fixed at a value equal to 
the characteristic measurement error. 

We determined the best-fitting model parameters and 
their uncertainties using a Markov Chain Monte Carlo 
(MCMC) analysis with a Gibbs Sampler (Geman & Ge- 
man 1984; Ford 2005; Winn et al. 2009). The MCMC fit- 
ting algorithm was implemented using the Transit Anal- 
ysis Package (TAP; Gazak & Johnson 2011, in prep.), a 
graphical user interface-driven analysis tool written in 
the Interactive Data Language (IDL). We constructed 
10 chains containing 10 6 links using initial conditions 
based on a simple least-squares fit to the phased pho- 
tometry. We chose step sizes such that 30-40% of the 
steps are accepted. We discarded the initial 10 5 links in 
each parameter chain to allow for "burn-in" before com- 
bining the chains. The resulting chains of parameters 
form the posterior probability distribution, from which 
we select the 15.9 and 84.1 percentile levels in the cu- 
mulative distributions as the "la" confidence limits. In 
most cases the posterior probability distributions were 
approximately Gaussian, and we therefore report only 
symmetric error bars for simplicity. The final, iterative 
fitting procedure and derived parameters are presented 
in § 6 

6. PHYSICAL PROPERTIES OF THE LHS 6343 
STELLAR SYSTEM 

We solve for the physical parameters of the LHS 6343 
system using the following iterative procedure: 

1. Fit the light curves to obtain the scaled semimajor 
axis, cir. We initially use AK P — 0.7, ok p = 0.1, 
AZ = 0.5, and oz = 0.1 in Eqn. 9 and refine these 
values with subsequent iterations. 

2. Estimate Mc based on our Keplerian fit to the RV 
time series by solving 



M c sin 3 i 



K 3 P 
2irG 



(M A + M C ) 2 V^ 



(10) 



We initially assume Mc <C Mb, and relax this as- 
sumption as Ma is revised. 

3. Use an and Mc to estimate Ma and Mb by mini- 
mizing Eqn. 8. 
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Fig. 4. — Comparison of the observed stellar and transit proper- 
ties (filled circles with error bars), and the corresponding quanities 
predicted by our model (red swaths; Eqns 2, 4 and 6). In order 
to show all parameters on the same scale, the quantities have been 
divided by the best-fitting model values. 
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Fig. 5. — Each panel shows the joint, posterior pdfs for two 
parameters at a time, marginalized over the remaining parame- 
ters. The contours show the iso- probability levels corresponding to 
{68.2, 95, 99.7}% confidence. 

4. Test for convergence in Ma, Mb, Mc and AKp 
by comparing the current values to those of the 
previous iteration. If the change is larger than 10% 
of the parameter uncertainties, then go to step 1. 

We find that convergence is rapid, requiring only three 
iterations. 

Our best-fitting stellar parameters are Ma = 0.370 ± 
0.009 M Q , M a = 0.30 ±0.01 M , d = 36.6 ±1.1 pc and 
[Fe/H] = 0.04 ± 0.08. The fit results in xLt = 2- 1 with 
8 data points and 4 free parameters, indicating that our 
fit is acceptable but that our photometric measurement 
uncertainties are likely overestimated. The best-fitting 
model values are compared to the observed quanities in 



Fig. 4. The marginalized, posterior probability density 
functions (pdf) are shown in Fig. 5. 

Our analysis of the light curve results in Rq/Ra = 
0.226±0.003 and a R = 45.3±0.6, with AK P = 0.74±0.10 
and AZ = 0.5 ± 0.1. The best-fitting transit model is 
shown in Fig. 1 for both the Kepler (blue) and Nickel 
data (red), along with the residuals. For both data sets 
we recover a r = 0, consistent with no red noise contam- 
ination. However, we do see additional structure in the 
in-transit residuals of the Kepler fit. We do not know 
the source of this increased scatter. Because there are 
so few points, and because the data do not appear to be 
time-correlated during any single transit event, we found 
that the extra scatter is accounted for by allowing the 
fitting procedure to increase the white-noise component, 
<j w . We find a w = 1.1 x 10 -3 (1.8-minute cadence) for the 
Nickel data and a w = 1.1 x 10~ 4 (29.4-minute cadence) 
for the Kepler data. 

6.1. The Radii of LHS 6343 A and LHS 6343 C 

Given the mass of Star A, we can estimate the stellar 
radius by evaluating Eqn. 7 which results in a stellar 
radius Ra — 0.378 ±0.008 Rq. The radius ratio from the 
light curve analysis, together with Ra yields the radius of 
LHS 6343 C, R c = 0.833±0.021 R Jup . The other physical 
properties of the brown dwarf are listed in Table 5. 

6.2. The Dependence of Rc on the Flux Contribution 
of LHS 6343 B 

In § 5.2 we describe our method of estimating the flux 
contribution of Star B, as parametrized by the magni- 
tude differences AA'p and AZ. These values are of cen- 
tral importance to the determination of the corrected 
values of a/RA and Rc / Ra as derived from the anal- 
ysis of the light curves, and hence the physical proper- 
ties of the brown dwarf. Unfortunately, we do not have 
spatially-resolved photometry of the binary stars, so we 
were forced to rely on theoretical stellar models to es- 
timate the magnitude differences. Owing to incomplete 
knowledge of molecular opacities, particularly in opti- 
cal bandpasses, the model grids provide only rough es- 
timates of the true flux contribution of Star B. To en- 
capsulate this imperfect knowledge in our analysis, we 
used relatively large uncertainties (0.1 mag) for AKp 
and AZ. These errors are propagated throughout our 
MCMC analysis and are reflected in the final physical 
properties of the LHS 6343 system. 

As a test, we performed independent fits to the Ke- 
pler and Z-band light curves. We found that the transit 
properties from the two analyses agreed extremely well, 
differing by only a fraction of a a in each value. Figure 6 
illustrates the dependence of Rc on AKp and AZ, and 
compares the value of Rc measured from the indepen- 
dent light curve fits. Based on this figure we feel that 
our value of AZ is reasonable since its value shouldn't 
be less than A J = 0.49. Similarly, and AZ shouldn't ex- 
ceed our estimate of AKp, and AKp shouldn't be much 
larger than our measured AV = 0.74. These results, 
along with the close agreement between the two mea- 
surements of Rc , provide additional confidence that our 
model-based magnitude differences are valid. 

6.3. Searching for Transit Timing Variations 




0.2 0.4 0.6 0.8 1.0 

A Z, A K p 



Fig. 6. — The dependence of Rc on AKp (upper, blue) and AZ 
(lower, red), based on the marginalized posterior pdfs. The con- 
tours represent the 68.2% and 95% confidence regions. The two 
filled circles with error bars show the final values of Rq, AKp and 
AZ from our analysis, and the error bars represent the correspond- 
ing 68.2% confidence bounds in each dimension. 

To measure the individual transit mid-times we fixed 
all of the global parameters (Rc/Ra, ci/Ra, P, % and 
the limb-darkening coefficients) and fitted each transit 
event separately using the MCMC algorithm described 
in § 5.3. Table 4 lists the time at the mid-point of each 
transit T m jd; the difference between the measured values 
and the those predicted by a linear ephemeris; and the 
formal measurement uncertainties which range from 10- 
23 seconds. We see no statistically significant timing 
variations. 

6.4. Limits on the Secondary Eclipse Depth 

To place an upper limit on the secondary eclipse depth 
we fitted a simplified transit model to the Kepler photom- 
etry near one-half phase away from the primary eclipse. 
We fixed the transit parameters from the fit to the pri- 
mary eclipse, along with e and U) the from the RV anal- 
ysis; assumed no limb darkening; and used final value of 
AKp = 0.74 ± 0.10. By allowing only the transit depth, 
a w and the terms describing the out-of-transit normal- 
ization to vary, we place a 95% upper limit on the eclipse 
depth of d < 6.5 x 10" 5 . 

6.5. Limits on the System Age 

Measuring the ages of M-type stars is notoriously dif- 
ficult, except for the rare cases when stars are in clusters 
or associations, or show indications of extreme youth. 
The space motion of LHS 6343 rules out membership in 
any known moving group or cluster. We therefore must 
rely on spectroscopic and kinematic indicators that can 
at least tell us if a star is either younger or older than a 
few Gyr. 

Our HIRES spectra show remarkably low chromo- 
spheric activity in the CallK line (Fig. 7). None of the 
140 low-mass stars on the CPS M dwarf survey have a 
chromospheric S values, as measured on the Mt. Wilson 
scale (Wright et al. 2004), as low as what we measure for 
LHS 6343 A. After attempting to correct for the « 30% 
contamination from star B to the emission core, we es- 
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Fig. 7.— HIRES spectra near the Call K line for LHS6343 and 
several M dwarfs with ages greater than 1 Gyr. 

timate S « 0.4 for LHS 6343 A. The most chromospher- 
ically quite M dwarf in the CPS sample is G1445, with 
an average of S — 0.5. 

LHS 6343 was not detected by ROSAT 12 , but at 36.6 pc 
this only rules out extremely young ages less than ~ 
100 Myr. Our HIRES spectra exhibit narrow lines in- 
dicating V ro t sin i < 2 km s _1 . No rotational modulation 
of star spots is seen in the Kepler light curve, indicating 
either a long rotation period (P rot > 40 days) or ex- 
tremely low spot coverage, or both. Thus, the age of the 
LHS 6343 system is most likely greater than 1-2 Gyr. 

7. DISCUSSION 

LHS 6343, C is remarkably similar to the brown dwarf 
NLTT41135B, recently discovered by the MEarth transit 
survey (Irwin et al. 2010). Both brown dwarfs orbit one 
component of a nearby M+M visual binary, both stellar 
systems are hierarchical triples, and both were discovered 
through transit photometry. However, owing to its more 
favorable orbital inclination, transit photometry provides 
a direct measurement of the radius LHS 6343 C. 

The large mass ratio between parent star and compan- 
ion argues against the formation of LHS 6343 C in the 
protoplanetary disk of LHS 6343 A. This is because there 
is likely not enough material in the disks of M dwarfs 
to build a 62.9 Mj up companion, whether through core 
accretion (Laughlin et al. 2004) or disk instability (Boss 

12 The Rontgensatellit (ROSAT) was a joint German, US, and 
British X-ray observatory operational from 1990 to 1999. 



T3 



1 1 1 1 1 1 1 1 1 1 


i i i | i 

OGLE-TR-123B 




1 








CoRoT-3b 


CoRoT-15b 










'****.. J ^"^-^^^ WASP- 


30b / OGLE-TR 


-122B " 


1Gyr 


_HS6343C / 




i ... i ... i 


i 




20 40 60 


80 


100 



Mass [M Jup ] 

Fig. 8. — The masses and radii of known low-mass, transiting 
objects. Also plotted are the predicted mass-radius curves from 
the Baraffe et al. (1998) interior structure models for 0.5, 1, 5 Gyr, 
top to bottom. Not shown are the young, eclipsing brown dwarfs 
in the 2MASS 2053-05 system (Stassun et al. 2006), which have 
radii well above the plot range. 



2006). It is therefore more likely that LHS 6343 C formed 
in a similar manner as LHS 6343 B, through the fragmen- 
tation of a portion of their natal molecular cloud. Irwin 
et al. (2010) make similar arguments regarding the origin 
of NLTT 41135 B. 

Prior to this year, the only other brown dwarfs with 
accurate radius measurements were the substellar com- 
ponents of the 2MASS J05352184-0546085 eclipsing bi- 
nary system Stassun et al. (2006), and Corot-3b (Deleuil 
et al. 2008). The 2MASS J05352184-0546085 system 
was discovered in a young star-forming region, and the 
two brown dwarfs have very large radii for their masses 
(0.669 R Q and 0.611 R ) because they are still undergo- 
ing gravitational collapse. 

In the latter portion of this year, two other transiting, 
field brown dwarfs have been discovered by transit sur- 
veys. The ground-based Wide Angle Survey for Plan- 
ets discovered WASP-30b (60.96 ± 0.89 M Jup , 0.89 ± 
0.21 Rj up ; Anderson et al. 2010), and the space-borne 
CoRoT mission discovered CoRoT-15b (63.3±1.1 M Jup , 

1.12±g;i S Rj up ; Bouchy et al. 2010). Both of these com- 
panions have masses comparable to LHS 6343 C, yet or- 
bit single, F-type stars that are much more massive 
than LHS 6343 A. The current sample of known transit- 
ing brown dwarfs are shown in Fig. 8, along with two 
very low-mass eclipsing M dwarfs. Also shown are the 
mass and radius predictions of the Baraffe et al. (1998) 
interior models for three different ages, illustrating that 
the radius of LHS 6343 C is consistent with the model 
predictions for a brown dwarf with an age > 1 Gyr yet 
< 5 Gyr. 

Our knowledge of the physical characteristics of brown 
dwarfs is starting to expand owing to the growing pro- 
ductivity and efficiency of wide-field transit surveys, both 
from the ground and in space. LHS 6343 is one of the 
closest stars in the Kepler field and points the way toward 
additional brown dwarf discoveries in the near future. 
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TABLE 5 

System Parameters of LHS 6343 



Parameter 


Value 


68.3% Confidence Comment 






Interval 



Transit Parameters 
Orbital Period, P [days] 
Radius Ratio, (Rc / R*,a)cott 
Transit Depth, (Rc /R*. a) corr 
Scaled scmimajor axis, a^ = a/R+ 
Orbit inclination, i [deg] 
Transit impact parameter, b 

Other Orbital Parameters 

Semimajor axis between Star A and C [AU] 

Eccentricity 

Argument of Periastron w [degrees] 
Velocity semiamplitude Ka [km s -1 ] 
Systemic Radial Velocity 7 [km s _1 ] 

Stellar Parameters 

M A [M Q ] 

M b [M e ] 

R*,a [Re] 

PA [p@] 

log 3a [cgs] 

[Fe/H] 

Distance [pc] 
Va [mag] 
V B [mag] 
(B - V) A 
(B - V) B 
AKp [mag] 
AZ [mag] 

Brown Dwarf Parameters 
M c [M Ju „] 

Rc [Rjup] 

Mean density, pc [g cm -3 ] 
loggc [cgs] 



12.71382 ±0.00004 A 

0.226 ±0.003 A 

0.051 ±0.001 A 

45.3 ±0.6 A 

89.50 ±0.05 A 

0.40 ±0.04 A 

0.0804 ±0.0006 B 

0.056 ±0.032 C 

-23 ±56 C 

9.6 ±0.3 C 

-46.0 ±0.2 C 

0.370 ±0.009 D 

0.30 ±0.01 D 

0.378 ±0.008 D 

6.6 ±0.4 A 

4.851 ±0.008 B 

0.04 ±0.08 D 

36.6 ±1.1 D 

13.88 ±0.03 E 

14.63 ±0.06 E 

1.57 ±0.07 F,G 

1.60 ±0.07 F,G 

0.74 ±0.10 F 

0.5 ±0.1 F,G 

62.9 ±2.3 B,C 

0.833 ±0.021 B 

109 ±8 B,C 

5.35 ±0.02 A 



Note. — Note. — (A) Determined from the parametric fit to the Kepler light curve. (B) 
Based on group A parameters supplemented by the photometric stellar mass determination 
described in § 5. (C) Based on our analysis of the Keck/HIRES RV measurements. (D) 
Based on our photometric mass and radius determinations described in § 5. (E) Eqn. 5. (F) 
Interpolation of Padova model grids. (G) (B — V) colors from the Padova models have been 
corrected by adding empirically measured offset of 0.21 ± 0.07 mag. Gunn-Z magnitudes 
were approximated using SDSS z'-band, with an additive correction of 0.07 mag based on 
the difference between the model A J compared to the measured A J in Table 1. 



